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Abstract 

Arguably, the greatest achievement of many-body physics in the fifties was that 
of developing the tools for a complete description and a thorough understanding 
of superconductivity in metals. At the basis of it one finds BCS theory and the 
Josephson effect. The first recognized the central role played by the appearance 
of a macroscopic coherent field usually viewed as a condensate of strongly over- 
lapping Cooper pairs, the quasiparticle vacuum. The second made it clear that 
a true gap is not essential for such a state of matter to exist, but rather a finite 
expectation value of the pair field. Consequently, the specific probe to study the 
superconducting state is Cooper pair tunneling. 

As forcefully stated by Phil Anderson, tunneling experiments not only gave a 
measurable physical reality to the phase of the pair wavefunction -the Josephson 
current is phase-dependent-, they also (essentially) eliminated any uncertainty 
concerning the mechanisms of electron-electron and electron-phonon interac- 
tions which are at the basis of the origin of pairing in metals. This also meant 
the end of superconductivity as a wide open, speculative field and the beginning 
of a thoroughly quantitative "exact" era, with uncertainties below the 10% level. 

From this vantage point of view it is not difficult to argue that important 
progress in the understanding of pairing in atomic nuclei -a subject of research 
started about the time of the BCS papers, but still far from having entered the 
"exact" era- may arise from the systematic study of two-particle transfer reac- 
tions starting from drip line, exotic, halo nuclei (like for example ^'Li and ^^Be), 
stabilized by the pairing correlations associated with a single Cooper pair, and 
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gradually proceeding to the case of many (but still few)-Cooper pair systems like 
e.g. Sn-isotopes. 

Time seems to be ripe for such a study, in keeping with the fact that one now 
knows how to describe microscopically the structure of these systems taking prop- 
erly into account the interplay between bare and medium polarization (induced) 
pairing interaction between Cooper pair partners. Also, how to calculate the ab- 
solute two-particle transfer cross sections taking properly into account the full 
non-locality of the Cooper pairs (correlation length much larger than nuclear di- 
mensions). While much of what is necessary to carry out the research program 
needed to enter the "exact" era of nuclear pairing making use of Cooper pair trans- 
fer is to a large extent known, a unified picture to implement it is still lacking. 

To remedy this situation and eventually help at starting the new rigorous era 
of pairing studies in nuclei we discuss in this paper both the past and the present 
of two-nucleon transfer reactions, as well as try to foresee the near future of such 
studies in nuclei, induced by light as well as by heavy ion reactions. 

Key words: pairing, finite many-body systems, two-nucleon transfer, tunneling 
PACS: 25.40.Hs, 25.70.Hi, 74.20.Fg, 74.50.+r 
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1. Introduction 



Two years ago physicists celebrated a ll around the world the 50th anniver- 
sary of BCS theory (|Bardeen et alJ (|l957ar ). lBardeen et alJ (|l957br) ). The success 
of such a theory is not so much, or better not only, the fact that it provided the 
definitive solution to one of the most spectacular phenomenon of all of physics - 
permanent (super)currents estimated to be stable for 10^" y- , but that it provided 
a paradigm for the pheno menon of s pontan eous symmetry breaking and associ- 
ated emergent properties (|Andersonl ([l994()). This paradigm demonstrated to be 
successful in a variety of fields starting from solid state physics and extending to 
nuclear and particle physics, field theory and astrophysics. 



Bohr, Mottelson and Pines developed, in the summer of 1957 (|Bohr et al 



(|1958|) ). the basis of the theory of nuclear superfluidity which was eventually ap- 
plied to the description of the nuclear stru cture, in particu lar to the calculation of 
the moment of inertia of deformed nuclei (iBelvaevI (119591)) an d of quadrupole vi 



brations of spherical and of def o rmed nuclei ( Bavman ( 1960a ). iBavman ( 1960b ) 
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In subsequent years the consequences of the phenomenon of spontaneous 
breaking of gauge symmetry in nuclei was investigated. In keeping with the fact 
that generalized rigidity in gauge space constitutes the most basic of the associ- 
ated emergent properties, two-particle transfer reactions are the specific probe to 
study t he individua l members of the associa t ed pairing rotational and vib rational 



bands (IBohd (Il964 . iBes and Broglial (Il966l) . iBohr and MottelsonI (119751) '). a feat 



which is totally out o f reach in the case of tunneling between two superconductors 



in condensed matter (|Josephsonl (|1962|) ). 

Due to the fact that the number of Cooper pairs participating in the nuclear 
condensate is small, one can study the phenomenon in terms of specific orbitals, 
some of which play an essential role in the transfer process (hot orbitals). For this 



must be measured as well as calculated ( 
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(1985 
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(|1991)). 



Much have been learned concerning pairing correlations in many-body sys- 
tems b y studying metallic grains a t low temperatures (c.f. OFarine and S chuck 
(l2002b. iFarine and Schuckl (1 1999b . IPerenboom eTaP d 198 lb . lAndersonl d 1959b . 
Lauritzen et al. (Il993)) as well as by using semic l assica l approximations (see 



e.g. iBengtsson and Schuckl (|1980l) . iKucharek et al.l (|1989l) . iBroglia and Winther 
(|2004|) and references therein). New perspectives in the study of pairing cor- 
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relations among nucleons are being o pened through the s t udy o f exotic nuclei 
lying along the drip lines in general (|Dobaczewski et alJ (119961) ). and of halo 
nuclei i n particular. Insight in t o thes e systems through t wo-particle transfer re- 
actions (ILenske and Schrieden (1 1998b. iKhan et al. (120041) ) , and also wit h high - 
energy knock-o ut reactions ( Hansen and TostevTrl ( 2003 ). Tostevin et al. ( 2004 ). 



TostevinI (|2007h ) is expected to shed light on the relative role the bare and the in- 



duced pairing interactions play i n regions of v ery low density (see e.g. Fig. 2 of 



Richted (|1993h . and Fig. 3.21 of lHeydd (| 19981 ). Br oglia, private communicat ion), 



but also on BEC of dilute Fermi gases (see e.g. iPethick and SmithI (120081) and 

references therein), where single-pair transfer experiments are not possible . 

The two-particl e transfer re actions carried out at G anil (|Chatterj iee etal.l(l2008b. 



Keelev et all (120071) ). TRIUMF (iTanihata et all (120081) ) and Dubna (iGolovkov et al. 



(120081) ) as well as those in the planning stage which eventually will be studied, 
are expect ed to expand the f rontiers of our knowledge concerning quasispin pair 
alignment (|Andersonl (|l958|) ) and dynamical pair correlations in nuclei to an ex- 
tent unthinkable only few years ago. 

In keeping with the insight pr ovided by the results of tunne ling experirnents 
concerning pairing in rnetals ( And erson ( 1969 ). |josephson (ll969h.lSca laDin ol (119691) . 
Mc Millan and Rowelll (fl969l) . lEsJ(ll973b . lOiaveil d 1973b . IjosephsonI (Il973b . 



Nambul ( 1991 )). much is expected to be learned concerning the effective, strongly 



renormalized interaction responsible for the presence of Cooper pairs in the nu- 
clear medium, by measuring the two-particle strength functions (energy depen- 
dence of the tunneling phenomenon). 

To be prepared to meet the challenge we remind, in the present paper, the 
basic experimental findings as well as the theory which is at the basis of the uni- 
fied description of nuclear structure and two-particle transfer reactions developed 
through the years, and draw the grand picture which is slowly emerging from the 
exploration of pairing in finite many-body systems, to the limit of a single Cooper 
pair. 



2. Basic experimental findings 

In this section we present some of the basic results obtained probing the atomic 
nucleus through two-particle transfer reactions. The information arising from 
both light and heavy ion reactions is briefly summarized. The bias towards light 
ion reactions reflects the experimental situation. 

Under the influence of Coulomb and of nuclear external fields which change 
angular momentum by two units (or by studying the y-decay of compound nuclei 
with high angular momentum content), nuclei respond collectively in terms of 
rotations and of vibration, as testified by spectra of the type shown in Figs. [Hand 
121 Not only the spectrum displays a simple behavior with angular momentum as 
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shown in these figures (quadratic in the first case, linear in the second), but also 
a large enhancement of the quadrupole matrix elements (as measured in terms of 
single particle units) connecting two successive states. 

In keeping with the fact that, as a rule, physical theories are based on the 
analogy whi c h one can es tablish between unknown things and simple phenomena 
(|Weill (|l966h . IWeill (|l999h ) is that the resuhs displayed in Figs. [3] and |4] are to be 
considered compared with those shown in Figs. [Hand [21 From this analogy it 
emerges that two-particle transfer data provide the experimental confirmation of 
the soundness of the co ncepts of pairing rotational and pai ring vibrational bands 
(|Bes and Broglial (|l966l) ). see also [Ring and Schuckl (Il980|) . Ch. 1 1, in particular 
p. 449. Measuring the associated cross sections in terms of pure two-particle 
configurations - equivalent to Weisskopf units in the case of Coulomb excitation 
and inelastic scattering- one obtains values of « 5 - 10 and « 25 - 50 in the cas e 



of pair i ng vibrational and pa iring rotational excitations respectively (iBroglia et al. 
(ll972l) . lBroglia et all (Il978h ) (cf. Figs. [3 [H [5] and [6]). 

While two-particle transfer reactions induced by light ions at energies of « 
5-10 MeV/nucleon, like {t,p) reactions, provide detailed spectroscopic infor- 
mation in terms of highly structured L-dependent angular distributions associated 
with the excitation of single nuclear states (see below), the two neutrons in the 
triton do not display the same correlations as a Cooper pair in nuclei. Among 
other things, this is t he origin of the over l ap factor Q„ ^ 1 app earing in the two 
particle form factor (jOlendenningl (|l968r) . iBroglia et al.l (119731 )). This limitation 
can be overcome by means of two-particle reactions induced by heavy ions. By 
appropriately choosing target and projectile one can prepare the Cooper pair to be 
transferred to have the desired properties, and tune the two-particle process to be 
the specific tool to probe nuclear pairing correlations. 

Arguably, the main drawback displayed by heavy ion reactions when used 
for nuclear structure studies is connected with the non-specificity of the angular 
distributions (bell shaped), associated with the very short wavelength of t he rel a- 
tive motion of the two ions (semiclassical behavior) (|Broglia and Wintherl (|2004l) ). 
Transfer essentially occurs at grazing conditions. For large relative distances the 
exponential decay of the Cooper pair wavefunction prevents transfer from occu- 
rring. At small distances many channels open up and the flux of the elastic and 
quasielastic channels is drained out (strong absorption). Thus the uncharacteristic 
shape of the angular distributions. 

The second major drawback of heavy ion reactions is the difficulty to observe 
single final states. The large Coulomb and hadronic fields induced by the mu- 
tual interaction of the two colliding ions makes transfer from and to inelastic and 
Coulomb excited states quite likely. 

In any case, experimental ingenuity combined with new technical develop- 
ments have made it possible to obtain precious information from these reactions. 



2 Basic experimental findings 



6 



30 



0) 20- 



152 

Decay of Dy 




30 40 
Spin (ti) 



60 



Figure 1 : A schematic representation of the proposed y-ray decay paths from a high-spin 
entry point in '^^Dy. The major initial decay flow occurs mainly via E2 transitions in 
the um-esolved y-ray continuum and reaches the oblate yrast structure between spins 30^ 
and AQh. A small 1 % branch feeds the superdeformed band, which is assumed to become 
yrast at spin 50?2- 55fi. The de-excitation of the superdeformed band around 26fi occurs 
when the band is 3-5 MeV above yrast, and a statistical type of decay takes into the oblate 
states betw een spin \9fi and 20?i The diagram also shows the low-deformation prolate 
band Cafter lNolan and TwinI dlQSSh ^ 
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Figure 2: Partial level scheme of ^'^CJ associated with the 



containing up to three phonons (After lAprahamian et al. 



uadru pole vibrational band 
1987h ). In the inset (up- 



per framed picture), a schematic representation of the energy spectrum of harmonic 
quadrupole vibrations as a function of the magnetic quantum number M is shown. 
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Figure 3: Pairing rotational band in even Sn-nuclei. Experimental energies of the 
J'^ - 0^ states of the even Sn isotopes. The heavy drawn Unes represent the values of the 
expression E = -B{S „) + Epyf- + 8.58N + 45.3 (M eV), where the binding energies B{A) (in 
MeV) are taken from lWapstra and Govd (|l97lh . The dashed line represents the parabola 
0.1 OCA'^ - 65.4)^. Also displayed is the excited pairing rotational band associated with the 
pairing vibrational mode. In all cases where more than one J" = 0"*" state has been excited 
below 3 MeV in two-neutron transfer processes, the energy cr(0,)£'(Ot)/ cr(Ot) of 
the centroid is quoted, as well as the corresponding cross section o"(0^). The quan- 
tity cr(Oj^) is the relative cross-section with respect to the ground state cross-section. The 
numbers along the abscissa are the ground-state (p, t) and (t, p) cross-sections normal- 
ized to the ^'^Sn <->'^^Sn(gs) cross-section. The (t, p) and (p , t) da t a utilizing in con- 



and (p , 

structi n g this figure were t a ken from refs. ( iBierregaard et all ( 19681) . . „ 
(Il969h. iFlvnn et all (Il970h . Ipleming et"aD dlQTOh l (After iBes and Brogha (1977h . and 
Brink and Broglial fcOOSh ). 



Bjerre gaa rd et al. 
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Figure 4: Pairing vibrations in even Pb nuclei. The ground state binding energies B(A) 
were used to evaluate the energies £'(A)exp as defined in the figure. The constant in front of 
the linear term, was determined so that £'(210)exp = £'(206)exp. Excited states were placed 
at the measured excitation energy above the relevant ground state. Experimental levels 
are shown as full lines, states predicted from the harmonic vibrational model ai^e shown 
as broken lines at an energy E(A)ih, whereas the short lines indicate the level position 
including a phenomenological anharmonic term. The two-particle transfer cross sections 
are given in terms of the pair-addition and pai r-substraction cross sections denoted x and 
y respectively (After Brink and Broglial ( 2005 )') 



In particular, based on the experience with light projectiles, an enhancement factor 
(EF) due to pa iring correlations for two-nucleon transfer betwe en heavy nuclei 
can be defined (|von Oertzen and Vitturil (|200lh . IWu et al.l (|l990|) ). This quantity 
is particularly simple to be extracted from the experimental data, provided the 
single-nucleon transfer probability Pi„ is used as the appropriate unit: P2„ = 
EFP\^ (see Figs. [51 |6l). In the independent particle picture this factor is just 
unity. Values of EF » 1 imply corr elations among pairs of nucleons transferred 
between the two heavy ions (see also lBroglia et al.l (|l972h ). 

Using particle-y coincidence techniques, well defined binary channels with 
final states of known quantum numbers can be selected. One- and two-nucleon 
trans fer cross sections in the collisio n of ^"^Pb-t-^'^Sn have been measured in this 
way (|von Oertzen and Vitturil (|200l|) ). Large enhancement factors have been ob- 
served, su ggesting a stro r ig con t ribution arising from sup e rfiuid pair transfer, a s 
predicte d dPietrich et all (Il97lh . IPietrich and Haral Jl973b . iBroglia et all (Il978b . 



see also Lotti et al.l 



mil)). 



3. The Josephson effect 

It is well-known that the phase of the wavefunction and the number of nu- 
cleons are conjugated variables. Gauge invariance, that is, invariance under phase 
changes, implies particle number conservation in the same way in which rotational 
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Figure 5: Measured transfer probabilities (transfer function) for one-(to;7) up to four- 
neutron (bottom) transfer in the system ''^Sn+'^°Sn, as a function of the overlap pai^am- 
eter do =Rminl{^]l^ + '^Y^)' ^mm being extracted in a unique way from the reaction angle 
for reactions at energies below or close to the Coul omb barrier Er. D ata at three energies 
(4.25, 4.55 and 4.8 MeV/nucleon) are shown. After lWu et aD(ll990h . 
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Figure 6: One-neutron (above) and two-neutron (lower part) transfer probabilities in 
206p|j^ii8^jj collisions as a function of the overlap parameter do measured in particle-y 
coincidences, with the Crystal Ball in anticoincidence. Filled symbols refer to a bombard- 
ing energy of 5.14 MeV/u, and open symbols to 5.32 MeV/u. The nuclei are identified 
by their characteristic y-transitions. The squared value of the measured In-transition 
probabilities are drawn as broken lines in the lower part. The shift between the squared 
of the l«-transfer probability and the corresponding 2»-tra nsfer probability defines the 
enhancement factor EF (after lvon Qertzen and Vitturil (|200lh ). Values of the order of 10^ 
ai^e observed in keeping with the fact that both target and projectile can be viewed as 
superfluid nuclei, their respective gs — » gs transitions bei ng enhanced by a factor of the 
order 25-30 as mentioned in the text (iBroglia et al. 
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invariance implies angular momentum conservation. A simple demonstration is 
the following: take a wavefunction with a fixed number of particles, 

W = alal---al\0). (1) 

Multiply the field operator by a phase factor 

Y = ie''^a\)ie'^al) ■ ■ ■ (e'^4)|0) = e'^^*P. (2) 

Clearly -id/d(f> = m". Thus 

* - (3a) 

i — = (3b) 
dN ^ 



and 



[^,A^] = /. (4) 



Correspondingly, there is an uncertainty principle AN Ac/) = 1. In keeping with 
this fact, the BCS solution of the pairing Hamiltonian 

\BCS) = Y]{U, + Vy"^alalm, (5) 

v>0 

violates gauge invariance, and \BCS ) describes a deformed state which defines a 
privileged orientation in gauge space (see Fig. |7]). The magnitude of the deforma- 
tion is measured by e"^a(, = (BCS \P^\BCS ), where = e^"^ Xiv>o ^v^J the pair 
creation operator. Of notice that the \BCS ) state can be viewed as a Bose-Einstein 
condensate of stro ngly overlapping Cooper pairs Q . Josephson proposed in 1962 



(|Josephsonl (|1962l) ) that there should be a contribution to the current through an 
insulating barrier between two superconductors which would behave like direct 
tunneling of condensed pairs at the Fermi surface from one superconductor to 
another. 



Quoting from lAndersonI ( ~969h commenting on the original BCS papers; 



. . . The idea that Py = ala^, the pair field, has a finite mean value is there, although 
oddly enough the statement that Bose-Einstein condensation of the pairs is irrele- 
vant is also there. The reason is that at the time Bose-Einstein condensation was 
thought as, of necessity, a condensation of preexisting particles, not as the appear- 
ance of a macroscopic coherent field, as we should understand it . . . 
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Each of the superconductors can be thought as a rotor in gauge space (see 
FigE]), weakly coupled to each other. The rate at which quanta are exchanged 
between the two superconductors is given by 

where Hcoupi- P^P ~ ^2,01^2,02^^^ ^^^g ^ g [Brink and Broglial JlOOSb . Appendix 



Ni = (-N2) = kHcoup,N^] = -L-:::^ ^ sin(2(0i - ^2) + 2d), (6) 



coup\ 

L) is the coupling Hamiltonian, d a phase related to the insulating barrier, while 
the rotational frequency of the rotors corresponds to the chemical potential of the 
superconductors (0 = j^dH/dN = This means that if there is a difference 
in chemical potential between the two superconductors, which can be obtained 
by applying an external voltage, there will be an oscillating current running be- 
tween the two superconductors (DC- Josephson effect). In terms of the voltage 
differential Vi - V2, one can write 

Ni ~ sin(^(yi - V2)t + 2d). (7) 
n 

This shows that the frequency of the oscillating current is determined by the a- 
pplied voltage, the carriers having a charge 2e. 

One can thus subsume superconductivity (as well as superfluidity) under the 
general theory of spontaneous symmetry breaking phenomena if one lets the order 
parameter he S = e^"^ao (note that A = G6 ; see however next section) and fix it in 
phase ((f)) as well as in amplitude (ao). Because of the long-range-order (10^ A < 
^ = fivpllls. < 10^ - 10"^ A), only a very small external force is necessary to fix the 
order parameter, and a direct observation of the phenomenon (experiments) can 
be carried out. 

The importance of the Josephson effect, then, is that it provides an instrument 
which can act in a similar way as a clamp for a solid or a coercive field for a 
ferromagnet: it can pin down the (relative) gauge phase through correlated pair 
transfer, the specific probe of dynamic (pairing vibrations) an d static (pa iring ro- 
tations) spontaneous symmetry violation of gauge invariance ([Anderson (fl984) ). 



Within this context we take issue at the otherwise superb monograph Quan- 
tum Liquids of Tony Leggett (Oxford University Press, Oxford (2007)). In the 
introduction Professor Leggett writes 

... it is neither necessary nor desirable to introduce the idea of "spon- 
taneous broken U{\) symmetry", that is to consider (alleged) quantum 
superpositions of states containing different total numbers of parti- 
cles; rather, I take from the start the viewpoint first enunciated ex- 
plicitly by C.N. Yang, namely that one should simply think, in non- 
technical terms, about the behaviour of single particles, or pairs of 
particles, averaged over the behaviour of all others . . . 
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In fact, it is this non technical thinking which tells us that in most nuclei, nucle- 
ons feel the average effect of all other nucleons in terms of a potential displaying 
a finite average value of the quadrupole operator (Nils son potential). Rotational 
bands (see Fig. [B observed up to the limits of the centrifugal force which the nu- 
clear liquid drop can sustain (/ « 60Ji) before fissioning, are some of the clearest 
evidence of the physical insight provided by such a picture. J. P. Elliott correctly 
argued, based on the results of his S U3 mod el that he could get rotationa l bands 
without the need of a deformed rotating body (lEllioti(ll958ah.lEllioti (|l958bl) , Elliot 
1962h . lElliot and Harveyl(ll963b . lBroglia and Maquedal(ll96dV)/ This is in princi- 



ple co rrect, and applies also to the Inter acting Boson Model (|lachello and Arima 
(Il988h . see also 
izations (see e.g 



Bes and Brqglial (|1979l) ') as well as to full shell model diagonal- 



Rowd (|1970|) and references therein). 



The power of the spontaneous symmetry breaking picture (see also iLaughlin 



(|2006h ) is that it can be used to systematically predict the collective degrees of 
freedom of many-body systems, (in particular of finite many-body systems of 
which the atomic nucleus is a much studied example) as well as to individuate the 
specific experimental probes to study them. The associated steps to be taken to 
implement this strategy are: a) find the variety of mean (constrained) fields which 
violate each of the different symmetries of the nuclear Hamiltonian (rotational, 
specular (parity), Galilean,isospin, gauge, etc.) b) restore symmetry by averaging 
out the alignment -privileged orientation in the different (abstract) spaces- thus 
giving rise to the variety of nuclear collective motions and associated generalized 
rigidities. Examples of such an approach are given by: 1) quadrupole and oc- 
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tupole rotations and vibrations (degeneracy of positive and negative parity states), 



2) dipole and 3) Gamow-Teller giant resonances, 
brational) bands H, etc (see e.g. Anderson ( 1964 ). 



(ll958h . lBohrand MottelsonI (119751) . iBes and Broglial(ll966h ) 



) pairing r otatio n al (and vi 



AndersonI (|1972l) . lAnderson 



The ground state wavefunctions associated with the variety of mean fields de- 
termined in step a) display a finite value of the: 1) quadrupole- and octupole-, 
2) center of mass-, 3)isospin- and 4) pair transfer-operator. These operators 
identify the experiments which specifically excite the corresponding collective 
nuclear modes, e.g. quadrupole or octupole Coulomb excitation in case 1) and 
two-particle transfer in case 4). 

The ratio between the energy of the pairing rotational transitions and of the 
intrinsic excitations (two-quasiparticle states) is much smaller than one (see Fig. 
[3]), the ratio approaching zero as the phase space available to correlate the nucle- 
ons, and thus the number of particles approaches the thermodyn amic limit. In this 
connection, it is illuminating to quote part of the discussion in IWeinbergI (|1996h 
on spontaneously broken global symmetries, where he uses a chair as an example 
of a macroscopic system: 

spontaneous symmetry breaking actually occurs only for idealized 
systems that are infinitely large. The appearance of broken symmetry 
for a chair arises because it has a macroscopic moment of inertia 3, 
so that its ground state is part of a tower of rotationally excited states 
whose energies are separated by only tiny amounts, of the order of 
. This gives the state vector of a chair an exquisite sensitivity 
to external perturbations; even very weak external fields will shift the 
energy much more than the energy difference of these rotational le- 
vels. In consequence, any rotationally asymmetric external field will 
cause the ground state or any other state of the chair with definite an- 
gular momentum rapidly to develop components with other angular 



^In other words, spontaneous symmetry breaking (2nd order phase transitions) is always as- 
sociated with emergent properties, in particular generalized rigidity. In fact, if one subjects a 
quadrupole deformed spheroidal nucleus to a torque, by pushing one pole (e.g. through an in- 
elastic {p,p') experiment), the opposite pole will react immediately as if the system was a rigid 
body, although the strong force through which the nucleons interact does not propagate instanta- 
neously. The moment of inertia of the associated rotational band (see Fig. [T]) gives a measure of 
this rigidity. The same applies to the tilting of one side of a domain wall defining a spin domain in 
a paramagnetic metal below the Curie temperature. In the phenomenon of spontaneous symmetry 
breaking, certain elements play a central role. For example, high angular momentum (hot) orbitals 
in the case of deformed nuclei and magnetic impurities in the case of spin domains in paramagnetic 
metals. In fact, the polarization effects of hot orbitals in inducing nuclear deformation (as a rule 
intruder states) can be studied in model calculations by artificially changing the value of spin-orbit 
potential, or simply by mutating (substituting) artificially the hot orbital by a low /-orbital. 



4 Successive and simultaneous transfer 



16 



momentum quantum numbers. The states of the chair that are rela- 
tively stable with respect to small external perturbations are not those 
with definite angular momentum quantum numbers, but rather those 
with a definite orientation, in which the rotational symmetry of the 
underlying theory is broken. 

Weinberg's arguments are true for an atomic nucleus only when N ^ oo. Never- 
theless, when one observes a (pairing) rotational band one can talk about a priv- 
ileged direction in gauge space, which can be clamped down by th e transfer of 

Coope r pairs from another super f iuid nu cleus in a heavy ion collision (IBroglia and Winther 
(120041) . Ivon Oertzen and Vitturil (|200ll) ). Because the total energy of the pairing 
rotational band is to be calculated in the laboratory frame, one has to add to the 
BCS Hamiltonian (Routhian) a term linear in A^, equal to AN, and associated with 
the Coriolis term in gauge space (equivalent to e.g. a»/, in the case of rotation of 
quadrupole deformed nuclei based on the intrinsic Nilsson state) to obtain it. The 
resulting energy of the pairing rotational band band AN + h^N'^/2^ contains thus 
a linear term in N and can consequently be viewed as the Anderson- Goldstone- 
Namb u mode (massless particle) of the associated gauge sym metry (lAnderson 
(Il958h . iGoldstonel d 196 ih . iNambul (fl96oh . iBroglia et aP dlOOOh l 

Of notice is that the same co ncepts have been instrurn e ntal to shed light int o 
other fields of phy s ics (see e.g. iNambu and Jona-Lasinid (119611) . IWeisd (119941) . 
Klimt et all Jl990l) . IVogl et all (fl990l) . IVogl and Weisd (Il99lb ) and of na tural re- 
search at l arge. In particu lar, into protein evolution and protein folding (jBroglia 
(120021) and lBroglial (120051)) and as a re s ult, in the design of drugs which do not cre - 
ate resistance (ILo Cicero et al. (120071 ). iFerramosca et al.1 (|2008b . iBroglial (|2009al) . 
Broglia (l2009bL iBroglia et al.l (l2009l)). In this case the second order phase transi- 
tion ( Shakhnovich and Gutinl ( 1993 ). Ramanathan and Shakhnovich ( 1994 )) takes 
place under evolutionary pressure (lowering the temperature in sequence space) 
and results in the becoming of proteins taking place when the "hot" amino acids 
of the protein are set in place. The most conspicuous emergent property is folding 
which is at the basis of metabolism and thus life on earth. Mutation of amino 
acids on cold sites of a protein gives rise to families of proteins, the associated 
native states displaying very similar energies (stability). They can thus be viewed 
as members of a rotational band in sequence space. The 2nd order phase tran- 
sition and thus one of its associated emergent properties, na mely folding ca n be 
blocked by mutating hot amino acids (protein denaturation) (IBroglial (|2009b() . cf. 
also footnote 2 of present section). 



4. Successive and simultaneous transfer 



Because the correlation length ^ of Cooper pairs in nuclei is of the order of 
30-40 fm, that is much larger than nuclear dimension (5-7 fm), one expects that 
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successive transfer of nucleon pairs does not destroy the correlation existing be- 
tween the nucleons forming the Cooper pair and contribute similarly than simul- 
taneous transfer to the associated two-particle transfer process . Furthermore, 
because the single-particle bases used to describe the transferred nucleons in en- 
trance and exit channels are, as a rule, different, non-orthogonal contributions are 
also expected to be important in the calculation of the absolute two-particle trans- 
fer cross section. It is of notice that only in the case of the two center shell model 
basis, two-nucleon transfer is tantamount to successive transfer, aside from the 
(very small) simultaneous contribution arising from transfer through the pairing 
force. 

Within this context one can mention that because the nuclear pairing inter- 
action has matrix elements which are considerably smaller than the value of the 
single-particle pot ential U(r) at grazing distanc es, the pair transfer potential looks 
much like {U(r))^ (|Broglia and Wintherl (120041) ). Wit hin this context i t is of inter- 
est to ment ion the disc i ission between J ohn Bardeen ( BardeenI (|1962h ) and Brian 
Josephson (IJosephsonl (|1962h ) (see also ICohen et al.l (|1962h ) concerning the fact 
that in a Josephson junction the gap vanishes at the surface of the insulating diox- 
ide layer and its implication concerning the possibility of a supercurrent across 
the barrier. 

Let us remind that in Gorkov's formulation of superconductivity, 
F(x, x')(~ 2v UyVy(py(x)(py{x') rclatcd to 6) is the amplitude for two electrons, at x 
and x', to belong to a Cooper pair. The gap function A{x) is given by V(x)F(x, x), 
where V(x) is a local two-body interaction at the point x. In the insulating bar- 
rier V(x) is zero and thus A{x) is also zero. The crucial point, however, is that 
vanishing A does not imply vanishing F (expectation value of the pair transfer 
operator). In fact, F(x,x') can have large amplitudes for electrons separated by 
distances \x - x'\ up to the coherence length. Hence, for barriers thin compared to 
that length, two electrons on opposite sides of the barrier can still be correlated. 
Consequently, the pair current is substantial, even if the mechanism of transfer is 
through a single-particle mean field potential. 

While in the nuclear case one would not observe supercurrents, one can study 
the consequences of dynamic and static deformation in gauge space in terms of 
the tunneling (transfer) of single Cooper pairs. Particularly illustrative is the case 
of the exotic nucleus ^^Li, embodiment of a single Cooper pair system. 

If neutrons are progressively added to a normal, stable nucleus lying along 
the stability valley, the Pauli principle forces them into states of higher momen- 
tum. When the core becomes neutron-saturated, the nucleus expels most of the 
wavefunction of the last neutrons outside to form a halo, which because of its 
large size can have lower momentum. In the case of ^'Li, one of the best studied 
halo nucleus, the last two neutrons are very weakly bound (by « 300keV). Con- 
sequently, these neutrons need very little energy to move away from the nucleus. 
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There they can remain in the "stratospheric" orbits, spreading out and forming 
a tenuous halo. If one neutron is taken away from '^Li, a second neutron will 
come out immediately (in keeping with the fact that '"Li decays into ^Li as fast 
as it is formed), leaving behind the core of the system, the ordinary ^Li. This 
result testifies to the fact that pairing plays a central role in the stability of '^Li. 
In fact, essentially all of the experimental findings characterizing ''Li can be ex- 
plained, even quantitatively, in terms of the following scenario: the formation of 
one Cooper pair which is held together to the ^Li core essentially by the coupling 
to low-frequency ph onons corresponding to long- wavelength vibrati ons of the 
core and of the halo ( Barranco et al. ( 2001 ). see also Vinh Mau ( 1995I)). A simi- 
lar picture accounts also for the nuclear structure of '^Be (|Gori et alJ (|2004t) ). Of 
notice, is that the exten ded shell mod e l resu lt s of this halo nucle us are consistent 
with the above picture (jSagawa et al.1 (1 1993b. iNavin et alJ (|2000f) '). as well as the 



dipole recoil model developed by lEsbensen et alj (|1997|) (see also iBroglia et al. 
(l2002b ). 

In the case of the Sn-isotopes, typical superfluid nuclei lying along the stabi- 
lity valley, the condensate is made out of about 5-10 Cooper pairs correlated half 
through the bare nucleon-nucle on interaction, the rem aining glue being provided 
by medium polarization eff'ects ([Barranco et al.l (|2004h ). 

The above picture of halo and of superfluid nuclei can be experimentally 
tested, in particular through processes of Cooper pair transfer induced by nuclear 
reactions. A particular exciting possibility is that of making "Li interact with e.g. 
'^''Sn, and measure the associated two-nucleon transfer cross section. Less spec- 
tacular, but equally interesting and simpler to analyze, is to measure reactions be- 
tween "Li and protons to determine the p robability that a proto n picks-up the two 
halo neutrons and form a tritium (t =1^2 (|Tanihata et al.l (|2008h ). see below). The 
correlation between the halo neutrons influences this probability, since it brings 
the two neutron which have to combine with the proton closer, as compared to the 
situation in which the two neutrons are assumed to be uncorrected (see Figs. [8] 
and|9]respectively). Of not ice also the two-neutro n transfer experiments from ^He 
in the reaction ^He-l-^^Cu (|Chatterjee et al.l (|2008b ). By analyzing such reactions, 
one would be able to obtain a measure of pairing correlations, and of the associ- 
ated pairing gap in low-density drops of nuclear stuff (halo nuclei), in a similar 
way in which (t, p) and (p, t) reactions on Sn- and Pb-isotopes have allowed to 
microscopically test the diff"erent contributions to the two-particle transfer cross- 
sections in stable, "spherical" and "deformed" nuclei in gauge space. 

In particular, changing the energy of the proton beam (in a (p, t) reaction), 
that is by measuring the two-particle transfer strength function one can, in princi- 
ple, match the frequency associated with the collective modes exchanged between 
Cooper pair partners (medium polarization effects). A resonant effect on the two- 
particle transfer cross section is thus expected (as observed in tunneling experi- 
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Figure 8: Spatial structure of two-neutron Cooper pair. The modulus squared wave func- 
tion |1'(n,r2)P = Kn,r2|0"'")P describing the motion of the two halo neutrons around 
the ''Li core (normalized to unity and multiplied by \6n^r^r^) is displayed as a function 
of the cartesian coordinates X2 - ri cos(0i2) and i/2 = r2 sin(0i2) of particle 2, for fixed 
values of particle 1 (ri = 2.5,5,7.5 fm) represented in the right panels by a solid dot, 
while the core ^Li is shown as a black circle. The numbers appearing on the z-axis of the 
three dimensional plo ts displayed on the left side of the figure are in units of fm~^ (after 
Barranco et aP (l200lh ). 
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Figure 9: Spatial distribution of the pure two-particle configurations s^^^ (0) and Pj^j*-^-* 
as a function of the x- and iy- coordinates of particle 2, for a fixed value of the coordinate 
of particle 1 {r\-5 fm). (after iBarranco et al.l (120011) ). 
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ments between superconductors (see e.g. Figs 32 and 33 of IScalapind (|l969t )). 
Of notice, that such an effect although not directly observed, has already received 
important circumstantial evidence from transfer experiments in stable nuclei. In 
particular, through the observation by means of one particle transfer reactions (as 
well as inelastic scattering), of the septuplet of states (hgi2(n) ® 3~)j resulting 
from the coupling of the h<)/2 prot on (ground state of ^"^B i) and the octupole vi- 
bration of the ^"^Pb core, see e.g. Bortignon et al. ( 1977b and refs. therein. The 
208p|^^3pjg j-j209gj experiments populating this septuplet can be viewed as the pro- 
cess equivalent to photoemission of e.g. C^q fuUerene and providing the electron- 
phonon coupling constants gy used to calculate the induced electron-electron in- 
teraction ari s ing fr o m the coupling of t he h^. electron to A^ , and Hg modes (see 



GunnarssonI (ll997h . lGunnarssonl (12004 . iBroglia et all (l2004 l 



In the case of the reaction ^^*'Po(?, a)^°^Pb((/z9/2(7r) ® 3~)) one obtains direct 
evidence of the renormalization of the single-particle state hg/2(n) by the 3^ vi- 
bration of the Pb core (see Fig. [TOl(c')'). 

Before concluding this section, let us remind some very general facts which 
are at the basis of nuclear superfluidity. The origin of the BCS many-body solu- 
tion of superconductivity can be traced to Cooper's solution of the two-electron 
problem in which all the other electrons play a role only through the exclusion 
principle. The ground state of ^^°Pb and of ^^Li, in which two neutrons move 
around the core of ^^^Ph and of ^Li respectively, provide nuclear realizations of 
Cooper's model. 

Much is known concerning two-nucleon pairing correlations of |^'"Pb (gs)). 
In keeping with the fact that this state can be viewed as the pair addition mode of 
^"^Pb (see Figs. |4]and[10l), it has been studied in detail in terms of two-nucleon 
transfer reactions. In this mode the two nucleon are correlated over a distance ^ = 
fi.vF/2Econ-, where Ecorr plays the role of the pairing gap for open shell, superfluid, 
nuclei. In the case of ^'"Pb, Ecorr ~ 1-2 MeV. Thus ^ = 25 fm. Of course, if the 
two nucleons are subject to an external field (the central potential generated by 
e.g. the ^°^Pb core), they cannot move away from each other more than 14 fm (see 
Fig. [TUl) . in keeping with the fact that the radius of ^"*^Pb is « 7 fm. On the other 
hand, in a heavy ion reaction with e.g. impact parameter 17 fm (see Fig. [111(a) and 
(b)), the central single-particle potential acting on one of the two nucleons to be 
transferred is much stronger than typical values of the pairing field. It will thus be 
this potential responsible for the transfer of one partner of the Cooper pair at time 
ti and of the second one at time ?2 (see Fig. [111(d)). And this two-step process 
will take place without loose of (pairing) correlation between the two nucleons. 
In other words, the Cooper pair is equally well formed at t < ti and t > (where 
the relative distance between the two neutrons is always less than 15 fm), than at 
ti < t < t2 where this distance can be much larger 24) fm. A similar argument 
applies for the non-orthogonality contribution (see Fig. [Tn(c)). 
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Figure 10: (a) Schematic representation of a two-paiticle con^elated state (pair vibration: 
pair addition mode) This object (Cooper pair) is highly non-local, the fermions partic- 
ipating in it feeUng the presence of the partner within distances inversely proportional 
to their correlation energy (light grey area) a fact which, in the case of halo exotic nu- 
clei, leads to an increase of the effective nuclear radius of the systems by factors of about 
50%. This may not sound much, if it was not for the fact that nuclear matter is highly 
incompressible and that small changes in the nuclear radius may imply nuclear instabil- 
ity. The mean single-particle field (dark grey area) can be viewed as an external field 
confining each of the member of Cooper pairs individually, (b) Pairing vibrations can be 
renormalized through the particle-vibration coupling mechanism in terms of self-energy 
and vertex connections, (c) direct one-paiticle pick-up can excite a 2p-\h like state while 
(d) two-particle pick-up may lead to a collective \p-\h excited final state. 
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Figure 1 1 : Schematic representation of the transfer of two correlated nucleons in first-order (c) 
and in second-order perturbation theory (d). In the first case a nucleon is transferred through the 
single-particle potential Uia acting at time fi (dotted vertical line). The second nucleon follows 
suit (dashed curve) through the non-orthogonality of the wavefunctions belonging to different nu- 
clei and, at time f2, it can be considered as a nucleon belonging to the nucleus B. In the second 
case, one nucleon is acted upon at f = fi by the single-particle field generated by the core A {Uia, 
dotted vertical line), and upon the second one at f = f2. In all cases the correlation between the two 
nucleons is maintained throughout in keeping with the fact that ^ ^ 30 fm. The wavy lines in both 
(c) and (d) indicate the correlations existing between the members of the Cooper pair, entity which 
is represented by two-particle arrowed lines. The light grey area in (a) indicates the (virtual) ex- 
tension of the Cooper pair wavefunction, the darker one the actual extension of the single-particle 
wavefunctions in the external single-particle field Uit- When the ion a comes to the closest dis- 
tance of approach from the target ion A (see (b)), the Cooper pair wavefunction turns from virtual 
to real extending now over essentially a length equal to ^. Transfer will receive (important) contri- 
butions over distances of ^ , as schematically shown with the help of a small volume of the Cooper 
pair wavefunction (solid circle). Of course, the transfer integral i'/cLper^^^^^'ficooper'^ receives 
contributions also from small distances, namely from all the dark grey area. 
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If we now concentrate our attention on ^^Li, arguably the best embodiment of 
a single Cooper pair fermionic system, in the intermediate, virtual state, the two 
neutrons can be as far apart as 6-8 fm from each other, and their center of mass dis- 
plays a mean square radius {rl^Y'^ « 5 - 6 fm (cf. Figs. [8l[9]and[10l). This in spite 
of the fact that the systematics of nuclear radii {R = 1.2A'^^) adscribes a value of 
« 2.7 fm to A = 1 1 . How does the system manage to extend so far away from the 
core? The answer is by binding very weakly the last two neutrons (two-neutron 
separation energy 52,, ~ 300 keV). Of notice that a gedanken experiment in which 
52„ is set to 20-30 keV would give rise to a Cooper pair which extends over tens 
of fm. This is consistent with the fact that the phonons (vibrational modes) which, 
exchanged between the two neutrons bind them to the ^Li core, are vibrations of 
the same mean field produced by this glue. If one knocks one neutron out, the 
other will follow suit, even if it is very far away (tens of fermis) from the other. 
There is nothing strange in this exquisitely non-classical quantum behavior. This 
is one of the reasons which makes Cooper pairs to be one of the most fascinating 
quantal objects found in nature (if nothing else because it changes the statistics of 
the particles involved, from Fermi-Dirac to Bose-Einstein). In the nuclear case, 
these fragile objects can be studied singly, through tunneling experiments (trans- 
fer reactions). Also at energies below that of the Coulomb ba rrier where they als o 
give rise to strong enhancem e nts of the fusion cros s sections (iDasso et al. (Il983b. 



Jacobs and Smilanskvl(ll983h . lBroglia et al.l(ll981h . lBroglia et al.l(ll983h . lBroglia 



(11985|) ). Pairs of nucleon separated by distances much larger than nuclear radii 
help each other in the (Cooper) tunneling through classical forbidden regions. 

Recently the {p,t) reaction on '^L i has been measured with the help of in- 
verse kinematics (|Tanihata et al.l (|2008b ). The differential cross section have been 
determined for transitions to the ^Li ground state and to the first excited 2^ state. 
Analysis of data has been carried out with success making use of correlated three- 
body wavef unctions. 

While final channel inelastic scattering events are likely to play a role in the 
^'Li(/), ff\A(2^) reaction, a direct two-particle contribution is expected within the 
model lea ding to the Cooper pairwavefunction shown in Fig. [H](see Fig. [TUl(d) 
as well as iBrink and Broglial (|2005h . Ch. 11). Within this context, it will be ex- 
tremely illuminating to calculate the absolute two-particle transfer cross section 
des cribing the Cooper pair wavefunction in terms of Nuclear Field Theory (NFT, 



see 



Bortignon et al.l (|1977h and refs. the rein), taking all possib le renormalization 



F K 

effects into account (cf. Figs. [SlfTOland lBarranco et al.l (|200l[) ). and the reaction 
process, within the framework of a quantal calculation (including successive, si- 
multaneous and non-orthogonality contributions, see Sect. 17]), and compar e the 
results with both theory and experiment as re ported inlTanihata et al.l (|2008h . Of 
notice is that NFT studies carried out on ^^Be (|Gori et al. I (l2004l) ) provide, through 
the analysis of single-particle spectroscopic factors, circumstantial evidence of 
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the central role the particle-vibration coupling mechanism plays in the structure 
of halo nuclei. 

5. Two-particle transfer reactions and time-reversal response function 

The fact that at the basis of BCS theory was not the pairing gap as initially 
thought, but the appearance of a coherent particle field which is distinct from it, is 
at the heart of the majority of results about superconductors. This can be seen by 
rewriting the pair transfer operator in the form 

where r is the time-reversal operator and Ivyj) = rlv/j). Equation ([8]) implies that 
the response of a superconductor to a time-reversal invariant perturbation is dra- 
matically different from the response of the system to a time-reversal noninvariant 
(pair-breaking or magnetic) one. Furthermore ([8]) explains why magnetic impuri- 
ties reduce the transition temperature and non magnetic ones do not. 

In the nuclear case magnetic impurities correspond to the intruder states (high 
j-values) lying close to the Fermi level and comin g down in energy from the up per 



shell due to the spin-orbit eff"ect. Backbending (iStephens and SimonI (|1972l) ) is 
associated with the breaking of pairs built on these states (typically iu/2(v) and 

Let us study the case of a nucleus under the effect of an external noninvari- 
ant time-reversal magnetic field (rotation). The model Routhian (Hamiltonian 
referred to the intrinsic, body-fixed reference frame) used to describe such a situ- 
ation can be written as 

H'^ = Hl + H,. (9) 
The first term describes the average rotating potential, 

Hl = Hl-oji,, (10) 

with 

(Nilsson Hamiltonian) Q being a quantity prop ortional to the quadrupole moment 



of the system (see e.g. iBes and Broglial (|1979r) ). 



The Hamiltonian Hf^ is invariant under space reflection (parity) and under 
rotations through 180° about the x-axis, i.e., rotations induced by the operator 

= expi-inj,), (12) 

but not under time reversal, with exception of the case in which co = 0. Conse- 
quently, the eigenstates \K) of //^ can be labeled by the parity (n = +1) and the 
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signature (a = +1/2) quantum numbers (jBohr and MottelsonI (|l974) '). Let us for 
simplicity consider particles moving in a single 7-shell, e.g. a cranked /13/2 orbital 
displaying axial symmetry. The eigenstates of ^ can be written as 



\j) = \n,a= 1/2) = Y,G'^(aj)\K) 

K 

\j) = \n,a = -1/2) = J]G{ico)\k) 



(13) 
(14) 



where the states \K) and t\K) = \k) are linear combinations of the two-fold de- 
generate Nils son states. 

Of notice that the states \ j) and \j) are related by the operation of time reversal 
only at ct> = 0. The violation of this symmetry is measured by the deviation from 
1 of the pair matrix element (coherent field) 



Mjj, = {jj\P^O) = J]Giico)Hl 



(15) 



The square of these matrix elements, in the case of a single j-shell (Z13/2) model, 
are shown in Figure [121 • While at a» = \j) = T\j) = \j) (Cooper pairs), at 
CO = l.O MeV all orbitals have become aligned with the rotation, leading to a 
complete violation of time reversal invariance. The critical rotational separat- 



ing the two extreme 
of these issues c.f. 



regimes is about c orot ~ 0.5 MeV. For a detai l ed discussion 



Nikam and Ringl Jl987h. iNikam et all (Il986h . iNikam et al 



(| 19871) . IVigezzi et al.l (119881) . iBrink and Broglial (120051) . and references therein. 

Making use of the language of linear response theory, we can view as an 
external field which adds to the system a pair of particles in time reversal states 
(see Eq. ([8])). One can study such a process as a function of co as shown in Fig. 
[T3] where we display a schematic representation of a heavy ion collision in which 
a deformed superfluid (target) nucleus is set in strong rotation by the Coulomb 
field of the projectile before the two-particle transfer process takes place. If the 
symmetry axis of the target nucleus is oriented at 45° or 135° the projectile exerts 
a maximum torque on the system. Oriented it at 0° or 90° with respect to the beam 
direction, essentially no angular momentum would be transferred to the deformed 
nucleus. Assuming that in the first situation co « cocr at the distance of closest 
approach, while a> « in the second case, will lead to a well defined pattern in 
the d^cr{A ^ A -I- 2)/dIdQ. dou ble differential cross section eventually determined 
making use of y-coincidence (|Broglia and Wintherl (|2004l) ). For co ~ 0, only the 
yrast yrast transition will be observed, the ymst (dizziest) band being the loci 
of the states of lower energy for a given angular momentum. For co > cOcr, that is 
after the superfluid normal transition takes place, one can distinguish between 
particle and hole states (see FigfT4l). In this case, we are dealing with an extension 



5 Two-particle transfer reactions and time-reversal response function 27 




Figure 12: Contour plots for dilferent values of oj of the two dimensional 7 x7 pairing ma- 
trices M\ connecting all the negative signature /-states with those of positive signature 
j. The soUd circles represent the square of matrix elements larger than 0.2. Consequently, 
the dashed zones indicate the region where pairing correlations are about one-half of the 
value displayed by the diagonal configurations at 6l» = 0. 



of the experiments exciting pair vibration bands to the case of finite cx), that is, a 
combination of Figs. 1 1131 and Fig. IH 
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Figure 13: Snapshots of a trajectory of the reaction ^^Kr+^^^U. The target nucleus is 
started with an orientation of 45° of the symmetry axis wit h respect to the incoming bea m 
so as to maximize the trans fer of angul a r mom enta (see iBroglia and Wintheii (120041) ). 
After basso et all (|l984l ) and basso etaP (|l982l) . 
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Figure 14: Schematic representation of the two-nucleon transfer processes exciting the 
one- and two-phonon pairing modes in the A'^ system starting from the corresponding 
yrast state of the N -2 system. The yrast yrast transition corresponds to the distinction 
of a pair removal mode, while the yrast (two-phonon) state corresponds to the creation 
of a pair addition mode. 
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6. Testing the tools 

6.1. Tunneling processes 

In connection with the micros copic descrip t ion of the tunneling process asso 



ciated with the Josephson effect, (|Cohen et al.1 (|1962|) ) one can use the Hamilto 



man 



H = Hi + H2 + ^ Tuq{al^a^^ + a^^^a^k^) + h.c. (16) 

Hi and H2 are the separate Hamiltonians of the two superconductors on each side 
of the barrier; Tkq is the exponentially small tunneling matrix element from state k 
on one side to state q on the other, and the relationship of phases shown is required 
by time reversal symmetry. 

One of the many procedures for arriving at (fT6l) is to find sets of single-particle 
functions for each side of the barrier separately, in the absence of the potential of 
the other metal: then one eliminates the nonorthogonality effect by perturbation 
theory. 

It is however necessary to be careful in using BCS perturbation theory in this 
case. This is because in a single, continuous block of superconductor it is quite 
correct to carry out the calculations using perfectly defined values of the relative 
phases of the energy gaps associated with the different states k and q, since the 
superconducting condensation energy actually results entirely from the precise 
coherence of these phases. It is however generally assumed that the total phase 
of the sample as a whole is meaningless. It is, on the other hand, not necessarily 
meaningless to discuss the relative phase of the two blocks of superconductors 
separated by an insulating barrier sufficiently thin for tunneling to occur. Clearly, 
again,the total phase of the assembly as a whole is not physical, but the relative 
phases can be given a meaning when we observe that electrons can pass back and 
forth between the two blocks by jumping across the barrier. Such a process leads 
to the possibility of coherence between states in which the total number of elec- 
trons is differently partitioned between the two sides -just as the phase coherence 
within the single block means that the number of electrons is not fixed locally 
(remember that ANAcp < In) and, for instance, there is a coherence between the 
state with N/l electrons in one half of the block and N/l in the other, and that 
with {N/l) -I- 2 on one side and (A'^/2) - 2 on the other. The same is true when 
referring to the left and right hand side superconductors in a Josephson junction. 



6.2. Two-nucleon transfer reactions 



Let us see how the above considerations translate in the case of the transfer of 
pairs of nucleons in a nuclear reaction, for example in a heavy ion reaction. 
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A very accurate picture of heavy ion reactions, in particular a two-particle 
transfer reaction, 

a{=b + 2)+A^b + B{=A + 2), (17) 
described by the total hamiltonian 

H = TaA + Ha + Ha + VaA, ^^^^ 

= ThB + Hh + Hb + VhB, 

is obtained by assuming that the relative motion of the centers of mass of the two 
ions can be described, in both entrance a(= (a. A)) and exit /3(= (b, B)) channels, 
classically (jBroglia and Winthed (120041)). 

In (fTSl) TaA is the kinetic energy of relative motion. Ha and Ha the Hamilto- 
nians describing the intrinsic degrees of freedom of nuclei a and A respectively, 
while VaA is the effective interaction between the nucleons in a and the nucleons 
in A. Similar notation has been used to describe the corresponding ions in the exit 
channel. 

All the information concerning the process (fTTl) is obtained by solving the 
time-dependent Schrodinger equation 

m-^ = Hip, (19) 
ot 

with the initial condition that the nuclei a and A are in their ground states, and that 
the relative motion is described by a narrow wavepacket xifp - Rpit)) of rather 
well-defined impact parameter and velocity. Expanding ij/ on the channel wave- 
functions ij/f} = il/''^(^h)ilfn(^B)ew(iSf)), according to 



^ - (20) 

Upon inserting (l20l) in (fT9l) . one obtains the coupled equations 

ifi ^ dp{,m^Wp)j^^^e-'''^'"' = 2<<Af I V, - Uy{ry)\ip,)^^ay{t)e-'^^'l\ (21) 

which, together with the condition ay{-oo) = 5{y, a) allows to calculate ap and 
thus the transfer reaction cross section, proportional to \ap\^. The sub-index on 
the matrix elements indicate that in the integration over the degrees of freedom of 
the two nuclei, the average center-of-mass coordinate r^y = ^{r^ + Vy) should be 
identified with the average classical coordinate, that is. 



ipy=\{ip+^y). (22) 
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Of notice that 

{iPy\Vy-Uy\iPy)=Q, (23) 

defines tlie ion-ion potential Uy to be the expectation value of the interaction Vy 
in the y-channel. 

A characteristic feature of the coupled equations (|2T)) is the presence of the 
overlap ((/^^|(^^) on the left hand side. If ^ and [3 describe two channels of the same 
partition, e.g. ^ = jS', the prime indicating excited states of the nuclei b and B, the 
overlap matrix is diagonal, namely, 

g{R) = («A^|<A/?)« = S(j3',/3). (24) 

If ^ and describe different partitions, the overlap {il/^\if/p) is different from zero 
in the region where the densities of the two nuclei overlap. 

The coupled equations (|2TI) can be written in a more compact way by intro- 
ducing the adjoint channel wavefunction 

^iJf = ^%Vr> (25) 

r 

fulfilling the orthogonality relation 

{co^,il,p) = 6{^,(3). (26) 
Solving (|2TI) in second order perturbation theory one obtains 

ap{t) = (a^(0)(0) + («/3(0)(i) + («/3(0)(2), (27) 



where 



{ap{t)\o) = ma\ (28) 

(«M0)(1) = 4 r - Ua\^a)jiJ'''''-''"^'''^dt', (29) 



and 



(a;,(0)(2) = (^) XJ^^'^^ " ^rl'ArH.ao^"''''"''^''''^^' 

«y —CO 



(30) 



The state vectors \co) have to include the non-orthogonality effects between chan- 
nels jS, y and a. To second order one finds 



2 {ilfaWy)^^^{llfy\ilfp)i^^llJa 



(31) 

-I- 
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Out of the very many second-order processes, we shall here only discuss the 
second-order effects in the two nucleon transfer reaction A+a(= b+2 nucleons) 
B(= A + 2 nucleons) + b, where the intermediate channel y corresponds to the one- 
nucleon transfer channel F(= A + I nucleons) + /(= b + I nucleons). 

Inserting oj^ in (|29|) and in (l30l) one finds that the two-nucleon transfer ampli- 
tudes can be written as 

a(oo) = (a;j)(i) + ia/j)orth + {.ap),ucc (32) 

up to second order of perturbation theory. The different quantities appearing in 
this equation are 

1 r°° 

(«Ai) = ^ dt{il,p\V,-U,\iP,)^J''^^'''^^''^\ (33) 

U —CO 

%J —CO 



\ ^ J(Ep-Ey)t/n 



i(Ey-E„)f/Tr 



(34) 



and 



{aB)orth = - ^ y, I dt{ifrMy)s 



Xco 
dt{i///3\if/. 
CO 



ymc^) (35) 

The first term (a^)(i) describes the simultaneous transfer of two nucleons. The 
two-step successive transfer is described by {ap)succ, while (afj)orth is a second- 
order contribution arising from the non-orthogonality of the wavef unctions asso- 
ciated with the different channels considered. 

Of notice is that the main contribution to (1331) arises from the single particle 
potential in channel a and the overlap between the single-particle wave functions 
of a nucleon in a(= b + 2) and in B(= A + 2). The contribution of the two-body 
pairing interaction, also present in (|33l) leads to a very small contribution. 

In keeping with the above discussion, in the independent particle model, that 
is in the case in which there are no correlation between nucleons 2y \^y){^y\ = 1^ 
(l35l) cancels exactly (1331) . The transfer reaction is then described as a purely suc- 
cessive transfer by the amplitude (l34l) . as expected. In the opposite limit of very 
strong correlations between nucleons one expects that two-particle transfer oc- 
curs essentially as a simultaneous transfer process. This can be seen rewriting 
(|33l) - (|35l) in the post-prior representation. In this representation, at variance from 
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the prior-prior representation of (1341) . the non-orthogonality term gets absorbed 
in the successive term which now reads 



%J —CO 

The expression (l33l) remains identical because of the post-prior symmetry. In the 
case in which the two particles to be transferred have a very strong mutual inte- 
raction Vu, («)2 become very small. This is because in the intermediate state one 
has to break a pair, an event which becomes less and less likely as Vn increases. 
Because in (1351) the transfer potential does not contain V^, the effect mentioned 
above implies that (a^)(2) ^ as oo- Now, in actual situations ^ (V), 

the first picture (i.e. a(i) ~ -{a)o,th) applies, and two-particle transfer process can 
be essentially viewed as a successive transfer process. 



7. Results 



7.1. Heavy ion reactions 

We now turn our attention to the results obtained in the analysis of the heavy- 
ion reaction 2°**Pb(^^0, '**0)^°''Pb at 86 MeV. As the wavelength of the relative 
motion is relatively short at this energy (A - 0.8 fm), the semiclassical scheme 
sketched above is expected to be well suited. In Fig. [15] we display the results of 
the semiclassical calculation for this reaction, in comparison with the calculation 
carried out within the framework of a fully quantal theory, (second order DWBA 
formalism) and the experimental data. 

In second order DWBA we need to calculate the simultaneous (T^^^), succes- 
sive {T Slice) and non-orthogonal {T)^q) contributions to the tr ansition amplitude 



between the initial {jjY and final (_//) single-particle states (|Bayman and Chen 



(119821)), 



T^'\ji,jf) = 2Y, f drfFdnidrA2UV^'(rAucr,)l'J'(rA2,cr2)tytB^^ (37a) 
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T^cUi, 7/) = 2 E f ^r/F^^rMrfrA2[^^>(rAi, cri)T^>(rA2, (^2)]°* (37b) 



o-',o-^ 



TfoUi, if) = 2 Z f ^r/F^rMrfrA2[^^'^(rAi, cri)^^>(rA2, cT2)t (37c) 

i' fl^ (Ti rTT 



/f, MO" 10-2 



x;rl-/(r,B)f;(rM)[^^>(rA2,c^2)^^'(rM,c^i)]^ 
X J Jr;iJr;2[^^^(r;2,cr^)^^'(r;i,^r;)]^ 



In these expressions, the spatial and spin coordinates of the two transferred nucle- 
ons are explicitly referred to with the subscripts 1 and 2. The subscripts A and b 
indicate the core to which the position of each of the nucleons are referred to. The 
vectors r^A, ^hB and Yff are the relative motion coordinates in the initial, final and 
intermediate channels respectively (note difference with Appendix B, where the 
intermediate-channel coordinate is denoted r^c)- 

The single-particle wavefunctions ^^(r, cr) with total angular momentum j 
are generated using a Woods-Saxon potential (see Table[I]) 

= \7T-V (38) 

l+exp[l^] 

where 

i?o = roA'l\ (39) 

The depth Vq is adjusted for each level to reproduce the corresponding experi- 
mental energy. The positive-energy wavefunctions ;\f[,'^^(raA) and;k'^~'(rfoB) are the 
ingoing distorted wave in the initial channel and the outgoing distorted wave in 
the final channel respectively. They are continuum solutions of the Schrodinger 
equation associated with the corresponding elastic scattering optical potentials. In 
Table[2]we show the (Woods-Saxon) parameters for these potentials. 

The transition potential responsible for the transfer of the pair is, in the post 
representation, 

Vfi = VbB - Up, (40) 
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Figure 15: Differential cross section for the reaction ^°^Pb('^0, '^0)^"^Pb at an energy of 
86 MeV in the laboratory frame, expressed in fib. We present the result s of the quantal and 
semicl assical calculations, along with the experimental data, see also Bay man and Chen 
(Il982h . 
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0.65 





1.24 
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Table 1: Parameters used for the single 
208pi3(i6Q i8o)206p|, dfiavman and Chen 



^artic le potentials associated with the reaction 
il982l) ). 
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where VhB is the interaction between the nuclei B and b, and is, the optical 
potential in the final channel. We make the assumption that VbB can be decomposed 
into a term containing the interaction between the cores A and b and the potential 
describing the interaction between b and each of the transferred nucleons, namely 

VbB = VbA + Vbl + Vb2, (41) 

where Vbi and Vb2 is the same mean field potential we have used to define the 
single-particle wavefunctions of the neutrons in the nucleus a. The transition 
potential is 

Ve = VbA + Vbl + Vbl - Up. (42) 

Assuming that {l5\vbA\a) - {/3\Up\a) (i.e, assuming that the matrix element 
of the core-core interaction between the initial and final states is very similar to 
the matrix element of the real part of the optical potential), one obtains the final 
expression of the transfer potential in the post representation, 

Vfi ^ Vbl + Vbl- (43) 

This last approximation seems reasonable when dealing with heavy ion reactions 
in which there is no charge transfer, but more care have to be exerted when deal- 
ing with reactions in which light ions are involved. To obtain the total pair trans- 
fer amplitude a sum of the contributions associated with each mean field state j 
(weighted with the correspondent spectroscopic amplitude Bj) is to be carried out: 

Tint = J] ^If^h (^''^On jf) + T^,,{h, ;» - Tf^{h, jf)) . (44) 

JfJi 

In Table |2] we show the single-particle active configurations and the cor- 
responding two-nucleon spectroscopic amplitudes associated with each single- 
particle level for the reaction under discussion. 

7.2. Light ion reactions 

We now turn our attention to two-particle transfer reactions induced with light 
ions. In this case, only a quantal calculation is possible. The calcula tion will 



be bas ed on formulae similar to those reported in Eqs. (I37al37bl37cl) (jBayman 



(Il97lh ). In Fig. [Hthe resuks for the ^^Ca(t, p)^^Ca reaction at 10.1 MeV, cal- 



culat ed in the post representat ion are shown in comparison with the experimental 



data (|Bjerregaard et al.l (| 19671) ). Of notice that the choice between either repre- 
sentation may not be arbitrary in this case, since the approximate cancellation be- 
tween the optical potential and the core-core potential is not equally well satisfied 
in each of the two representations. According to this consideration, in (t,p) and 
(p, t) reactions the representation corresponding to the proton-nucleus channel 
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V 


rv 


ay 


W 


rw 


61 


incoming channel 


100 


1.26 


0.45 


65.4 


1.26 


0.45 


outgoing channel 


65 


1.35 


0.34 


45 


1.34 


0.33 



Table 2: Paramet ers of the optica l pot entials associated with the reaction 
208p^^(-i6Q i8Q)206p|, (|Bayman and ChenI (jigsA ). obtained by fitting elastic scattering 
data. 



Nuclei 


n, I, j 


B 


206p|^ 


2;? 1/2 


-0.75 




1/5/2 


-0.47 




2p3/2 


-0.28 




0/13/2 


-0.28 




1/7/2 


-0.18 




0/19/2 


-0.14 


180 


0^5/2 


0.89 




2p3/2 


0.45 



Ta ble 3: Single-particle states and B-coefRcients (two-particle spectroscopic amplitud es, 
Brogha et all (Il973b ^ for the reaction 208pb(i6o,i8O)206pb dMaghone et all dlQSSh v 



see 



seems to be better suited to satisfy the cancellation mentioned above. In keeping 
with this reasoning, in the case of the reaction '^'^Ca(t, pY^Ca, the post representa- 
tion should be the preferred one. 

The triton wave function used in this and the '^°Sn(p, O'^^Sn reaction can be 
written as (see Appendix A) 



^{rip, r2p, ru) = p{rip)p{r2p)p{ri2)b(:icr,uMo-n2)]oXm,{o-p)- 



(45) 



The corresponding central single-particle potentials used to calculate the two- 
particle transfer form factors are displayed in Fig. \T2\ 

In Fig. [T8]the results for the ^^~Sn{p, tY^^Sn reaction at 26 MeV, calculated in 
the prior representatio n are displayed in comparison with the experimental data 
(iGuazzoni et al.l (|2006h ). Following the argument exposed in the discussion of the 
'*°Ca(f, p)'^^Ca reaction, in this case we should prefer the prior representation. Of 
notice that in all the calculations of the absolute cross sections discussed above 
there are no free parameters adjusted to reproduce the data. 
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Figure 17: Central single-particle potenti als of Ca ( solid l ine) and t (dashed Une) used in 
the quantal calculations shown in Fig. [T6](lBaymanl (1197 ih ). 
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8. Conclusions 

In a metal, the repulsion among the electrons compensates, to a large extent, 
the attraction generated by the positively charged ions. Such a system does not 
display the right cohesive energy. However, electrons being fermions, are de- 
scribed by an antisymmetric wavefunction. The associated exchange Coulomb 
interaction is attractive and adds to the cohesive energy an important contribution. 

This system,being localized in space and displaying rigidity, violates trans- 
lational invariance. One can observe it by shining a beam of light at it. Two 
outcomes are possible: a) electrons and ions oscillate rigidly together (Thomson 
scattering, center of mass oscillation), b) the system absorbs the photon which 
sets electrons and ions in oscillation against each othe r (photoabsorptio n, plas- 



mon). The plasmon state (Higgs' boson of the system, lAndersonI (| 19941) ). being 
orthogonal to the center of mass motion state -thus displaying 3A^ - 3 degrees of 
of freedom, A'^ being the number of electrons of the system-, leads to restoration of 
translational (Galilean) invariance. The emission and absorption by the same elec- 
tron, and the exchange between different electrons of virtual plasmons, gives rise 
to dressed particles which repel each other through a strongly screened Coulomb 
potential. 

A crystal not only violates translational invariance globally (center of mass of 
the system), but also locally in each point of the crystal where an ion is localized. 
Lattice vibrations -so called phonons (Goldstone or, more precisely, Anderson- 
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Goldstone-Nambu modes)- lead to an overall (translational) symmetry restora- 
tion. The dressing of electrons by phonons -and plasmons- is central to account for 
the observed specific heat of metals, while the coupling of electrons and phonons 
is at the basis of electric resistivity. 

The exchange of virtual phonons between electrons moving close to the Fermi 
energy gives rise to a tiny attraction, sufficient to overwhelm the tiny repulsion 
associated with the screened Coulomb field. The resulting overall attraction gives 
rise to pairs of electrons -Cooper pairs- of spin zero. At room temperature (« 25 
meV), these (quasi) bosons, being bound by an energy 2A of only a few meV, 
break as soon as they are formed. Eliminating thermal fluctuations by lower- 
ing the temperature of the crystal to few degrees K, allows for: a) condensation 
of strongly overlapping quasi bosons (Cooper pairs), b) distortion of the Fermi 
surface over an energy range 2A. This last effect selects a privileged orienta- 
tion in gauge space and thus leads to a spontaneous breaking of gauge invariance 
(non-conservation of particle number). Symmetry restoration is obtained by fluc- 
tuations of the "deformed" nucleus orientation in gauge space. In fact, pairing 
rotational bands are the fingerprints of deformation in gauge space. These bands 
have been observed in the case of finite systems (atomic nuclei). However, in the 
case of a metallic superconductor, the moment of inertia associated with pairing 
rotational bands is so large that it is not possible to observe the individual states of 
these bands. On the other hand, setting two superconductors in weak contact with 
each other, with the help of a thin dioxide layer (through which electrons can tun- 
nel) one observes, upon closing the circuit -that is, connecting among themselves 
the far edges of the superconductors opposite to the layer- a supercurrent (at zero 
bias). The carrier of this supercurrent has charge 2e (DC Josephson effect). This 
is a direct consequence of the coupling, through Cooper pair transfer, of the rotors 
in gauge space, associated with each of the two superconductors. In other words. 
Cooper pairs transfer gives rise to a coupling between the two superconductors. 
The detailed study of such processes have provided a complete picture of pairing 
correlation in metallic superconductors. 

While in condensed matter the actual calculation of the intensity of the Joseph- 
son current meets the formidable difficulties of a microscopic description of the 
tunneling across a dioxide layer, absolute two-particle transfer cross sections in- 
duced in nuclear reactions can be calculated with good accuracy obtaining an 
overall account of the experimental findings. The study of two-particle transfer 
nuclear reactions not only is expected to be instrumental in shedding light into 
pairing correlations in atomic nuclei, but also to gain detailed information con- 
cerning the relative role successive, simultaneous and non-orthogonality transfer 
play in the transfer process, as well as the role hot orbitals have in such processes. 
Also, to provide information concerning pair correlation across the barrier as a 
function of the di-neutron separation energy, from tens of MeV in stable nuclei to 
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hundreds of keV in halo nuclei. 
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A. Simultaneous transfer 



The ap pendixes A and B re l ate to light ion reactions and follow closely lBayman 



(Il97lh and lBavman and ChenI (Il982h . 

The simultaneous contribution to the (t,p) transfer cross section is 

in,mp 

where a summation is made over the spin orientations mt,mp of the incoming 
triton and the outgoing proton respectively, and the factor 1 /2 accounts for the 
two possible values of this orientation. The transfer potential Vir^i) is the proton- 
neutron potential acting on neutron 1, so the above expression is written in the 
post representation. 

For a A(t, p)B reaction, the relative motion between the nucleus A and the 
triton is represented by an incoming wave. If the two transferred neutrons are in 
the 5=0 singlet state, the spin is entirely due to the spin of the proton. We can 
thus write the incoming distorted wave as 



V / f\ / . y/47T(2l, + 1) 
rm'i^Aa, k,-, (Tp) = > exp (/crj gijXrAa)Yf;irAa) r Xm,(o-p)' (47) 

h 

where r^c is the vector between A and the center of mass of the triton, k, is the 
incident momentum, o-p is the spin coordinate of the proton and is the orbital 
angular momentum of the relative motion between the triton and the nucleus A. 

We have used Y'^ikj) = i'' -^^|^5m,,0' as k,- is oriented along the z-axis. Note that 
we have used the time-reversal phase convention, where the spherical harmonics 
are defined in terms of the associated Legendre polynomials as follows: 



V 4;r (/ + m) ! 

If we write 

Yl;(tAa)XmXo-p) = J]{l,0\/2 m,\jr m,) [Y''(tAa)x((rp)]^^ , (49) 
jt 

we can express the incoming distorted wave in the form 



,(+), , X V /• t\ V4;r(2Z,+ 1) 

i„j, "'^^^ (50) 



X 



(l, 1/2 m,\j,m,)\Y''(rAaMcrp)Y' . 

L ' Jm, 
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The outgoing proton distorted wave is 

(51) 

where Vgh is the position of the proton with respect to the nucleus B and cTp is its 
spin coordinate. Now, 

^ Yl:,{rBb)Yl::{Vif)x,n^{cTp) = J] Y',t(^f){lp m 1/2 m„\jp m + irip) 

X \Y'''(rBh)XmMp'^Y" (52) 
= J] Y',C,„^X^f){lp m - 1/2 m^Up m) [Y'^{fBb)Xnr„{cTp)]^ , 

and, finally, 

^|-;(r,,,k^,(r,) = ^Yj ^■'"e^P(-0/w/'-«^)^t"v(k,) 

f l„ip,m (53) 

Xilpin- nip 1/2 mp\jp m) [Y^''{rBb)xi(^p)Xl ■ 
We now turn our attention to the evaluation of the matrix element 

X (Ip m- nip 1/2 mp\jp m){lt 1/2 mt\it nit) i"'" exp[?XcrJ^ + o-'i)] 

nO* 



2Ytm,A) y r ^^^^W/,;v(rAl)W/;,>(r^2)[l^'^(fAl)l^'^(?A2)]!!^ 

><fipU{rBt)gi,j,{rAa)[x{(r,)x{cr2)^* [Y'''{fBh)x{crp)\^* V{r,,) 



with 



m, 

(54) 



r -^A2 - Tai 

1 

s =2 (fAi + rA2) - 

1 (55) 
n =2 (^Al + rA2) 

_ r^i + rA2 
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where r^i, rA2, r^ are the positions of the two neutrons and the proton with respect 
to nucleus A and cri,o-2 the corresponding spin coordinates. The wavefunction 



u,,j,{rAi)u,^j,(.rA2) [Y'firAi)Y'f(rA2)]l* (56) 

describes the two transferred neutrons in a final state with orbital angular momen- 
tum If and total angular momentum jf, and ^^(r, s) is the spatial wavefunction of 
the triton. 

The sum over 0-1,0-2 in (l54l) is readily found to be 1. We will now simplify 
the term [}''''(fB/,)A^(crp)j^'' [r''(rAa)^(o"p)]^' , first noting that, with time-reversal 
phases, 

[YHrBbMcrp)fj = (-l)i/2-P+i.- [YHrBbM-crp)]i^^ . (57) 
On the other hand, 



JM 

X ^YHrBb)x{-crp)]" [Y^rMcTp)]'^ ^ 



(58) 



In order to survive the integration, the angular and spin functions must couple 
to L = 0, S = 0, 7 = 0, so the only term that remains is 



{jp - m j, m,|0 0) >^Y'''{rsb)x{-(rp)\'' [Y\rAa)x{crp)Y^^ 
(-!)■ 



. ,;^0 (59) 



We couple separately the spin and spatial functions: 

{[y'(fB,);t(-(rp)]' [y'(fA«);r(t^p)]'}[ 

= ((/i)//i)^.|(//)o(ii)o)o \xi.-^p)xi^p)\ ^\^Bb)Y\rf,,^^ . 



(60) 
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We substitute (I57]),<l59l),(l60l) in ^ to obtain 



X (l nit - nip 1 /2 niplj m,){l 1/2 m,|j m,) r' exp[?'(crj' + crj)] 
X fij{rBh)gij{rAa) [Y\rsb)Y\r^a)'\^ y{n,)(P^{r, s) 

(61) 

The last sum over cTp is 

<T,, CTpin 

XXm(-0-p)X-ni(0-p) (62) 

The 9 j symbols can be evaluated to find 



((/i)//i),|(//)o(H)o)o = 



2j/ + 1 

2(2Zy + 1) 

27+1 
2(2/+ 1)' 



so 



(4^y 2j>+l 



(63) 



X </ - iTip 1 /2 mpl J m,){l 1 /2 | j ) z exp[/(crf + crj)] ^^^^ 
X V2y,^,,_„,^/k^) J ^^£^i,,^.^(r^i)i^/,i,(r^2)[i''^(rAi)>''^(fA2)]°* 
X MrBh)gij(rAa) [Y'(tBh)Y'(rj,,)]l Vir^Mr, s). 
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We now check the possible values of the Clebsh-Gordan coefficients, finding, for 
7 = Z-l/2: 



</ mt-nipl 12 mp\l-\ /2 m,)</ 1 /2 m,|/ - 1 /2 m,) 

I 



21+1 



if nit = m 
if m, = —m 



27TT ' " " (65) 



and, for j = I + 1/2: 



(l nit - mp 1/2 mp\l + 1/2 m,)</ 1/2 m^lZ + 1/2 m^) 

/+ 1 



— ifm. = m, (66) 



V?(m) 



Substituting, we get 



X 
X 



X V2y:,.,„^/k^) J ^^£^M,^.^,(r^i)M/,,>(r^2)[l^'^(?Ai)l^'^(?A2)]||* 

y(rM)0;i(r,s)[y'(fs,)y'(fA„)]|| 

[(///+l/2('"Bfo)e'//+l/2('"Ao)(^ + 1) + fn~\l2i.rBb)gU~\l2{rAc)l)6m,„m, 
+{fll+l/2(fBb)gil+l/2irAa) V^CTTT) - flI-l/2(rBh)gn-l/2(fAa) + l))^,,,^,-,,,,] • 

(67) 

We can further simplify this expression using 

[F'/(f^i)F'/(f^2)]°* = [Y'f{rAx)Y'f{rA2)^^ = J]{lf m If - m\0 0)7-{(f^i)Fi„(fA2) 

m 

= J]i-l)'^~'"{lf m If - m\0 0)Yl{(rAi)Y',f:(rA2) 

m 

= / yY!,UrAr)Y'S{tA2) 
V(2// + 1) 

(68) 
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where 6a is the angle between r^i and Tai- Now, 

[Y'(rBt)Y'(rAa)]l = J](iml - m\0 0)Yl(rBb)YUrAa) 

1 ^ (69) 
= -== y(-l)'^'"};UrB.)yi„,(f,J. 
V(2Z+ 1) ^ 

It can be easily shown that the integral of each one of terms in the above sum is 
independent of m, so we can drop the sum and multiply by 2/ + 1 the m = term, 
leaving 

[Y'(rst)Y'(^Aa)]l ^ (-1)' V(2/+l) Y',(rsb)Yl,irAa) ^^^^ 
= V(2/ + l)Y',(rBb)Y'*(rAa). 

We now change the integration variables from {VBb, r, »?) to (tau, a,P, y, rn, rip, r2p) 
(see Fig. [IS), 

d(r, T], VBh) 



ruripr2psmfi (71) 



d(rAa,a,fi, y,ri2,rip,r2p) 
being the Jacobian of the transformation. Finally, 

{¥-\kf)\v(rbi)\¥;\k,z)) = ^Yj V^Tr^^pt'K + 

, ^ r , r da dBdydr 12 dr I r,dr2r, sin B , 
X i;!„-„,,(k,) Jr^,y^(f^.) ^ ' ^ ^ -Y'oihb) 

X uifjf{rA])uifjf{rA2)V{rM)^Q{r, s)/'/,.(cos 0A)ri2r\pr2p 

X [(///+l/2(''efo)e'//+l/2(''Afl)(/+ + fll-\/2(rBb)gil-l/2(rAa)iym,„m, 

+{ftl+l/2(fBb)gU+l/2{rAa) ^Jlil + 1) - fll-l /2(rBb)gil-l/2(rAa) ^Kl + l))^mp,-m,] • 

(72) 

We note that the inner integral is a function of r^„ alone, and that it transforms 
as yo(f^f,) under rotations, because all the dependence on the orientation of Yac, is 
contained in the term Y^{rBb) ■ The inner integral can thus be cast into the form 

A(r,jy!,ir.j = f,am,,,r,,,r,,.in0 ^^^^ 

X F{a,/3, y, r^, rip, r2p, Rx, Ry, Rz)- 
To evaluate A(rAa), we put r^f, on the z-axis 

A{rAa) = 2nr' ^ J ^13 dy dri2dripdr2p sin^S ^^^^ 

X F(a,y8, y, rn, rip, r2p, 0, 0, r^J, 
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where a factor In has been included as the resuh of the integral over a, since the 
integrand clearly does not depend on a. We substitute (1731) and (I74l) in (I72l) . and 
after integrating over the angular variables of r^o, we obtain 



X ^i,-m„(k/) J"i^rA„ JjS Jri2 Jrip dr2p Tac, sinyS runpr2p 
X M;^y/rAi)M/,.y/rA2)V(r/,i)0o(r, s)i'//cos 0a)^/(cos %) 



X 

+( 



[{fn+i/2(rBb)gii+U2(rAa)(l + 1) + fii-\i2irBb)gn-\i2irAa)l)5m,„m, 

-{fn+\i2{rBh)gii+\i2{rAa) 1) - fn-\i2irBb)gii-i/2(rAa) l))^,,,,, /''sfe, 

(75) 

where we have used 



2/+ 1 
4;r 



-PKcos^). 



(76) 



The final expression of the differential cross section involves a sum over the 
spin orientations (see (|46|)). When rtip = 1 /2, m, = 1/2 or = -1/2, = -1/2, 
the terms proportional to 5mp,m, will include the factor 



when rrip = -1/2, nit = 1/2 

|y,t,_„,/ky)5.^,_„j = |y|(k^)i = 

and when nip = 1/2, = -1/2 



21 +\ 
An 



(77) 



I2/+1 1 
An 1(1+1) 



(78) 



12/ +1 1 
An 1(1+1) 



Pjicos 6) 



We use the tritium wavefunction 

0o(r, s) = p(rip)p(r2p)pirn). 



(79) 



(80) 



p(r) being a Tang-Hemdon wave function (ITang and HerndonI (Il965|) ). We obtain 

1 



• f 



'^(iC(«)p+iC(«p). 



(81) 



A Simultaneous transfer 



50 



with 

I 

X ^ dvAa d/3 dydrn dripdr2p ^.mp p{rip)p{r2p)p{rn) (§2) 

X Uifjf{rAi)uifjf{rA2)V{rb\)Pif{co% 9a)Pi(.cos 9if)rnripr2p 

X [fn+\i2irBb)gii+\i2irAa) (/+!) + fii-i/2(rBb)gii-i/2{rAa) i)lrBb, 



and 



X J" drAad/3dydri2 drip dr2prAa sin fip{rip)p(r2p)piri2) (-33) 
X M/^;y(rAi)M/^y/rA2)V(ri,i)i'/^(cos 0a)^/(cos 6if)ruripr2p 
X {flMI2{rBb)giMI2{rAa) - fn~i/2(rBb)gii-i/2{rAa)) / rsh- 

Note t he absence of the (-l)j factor with respect to what can be found in Ib ay man 
(|l97lb . due to the use of time-reversed phases instead of Condon-Shortley. This 
is compensated in the total result with the same difference in the expression of 
the spectroscopic factors. This ensures that, in either case, the contribution of all 
the single particle transitions tend to have the same phase for superfluid nuclei, 
adding coherently to enhance the transfer cross section. 

The above expressions can easily be rewritten for a heavy-ion reaction be- 
tween spinless nuclei, the distorted waves in entrance and exit channels being 



iff^^\rAa, ^Aa) = V exp (ia]) ^,4(r„A)^^^^^^^, (84) 

and 

<A<-)(r,B, k;,B) = ^'^ (-^'^f) f-S^^^^ Z '^l^^bB)Y^nihB). (85) 

If we are dealing with a heavy ion reaction, ^^(r, s) must be replaced with the 
spatial part of the initial two-neutron wavefunction 

»F(rH, rfe2, (Tx,(T2) = [iA^''(rM, cri)i/r^''(rfe2, 0-2)]° , (86) 

where r^i, r;,2 are the positions of the two neutrons with respect to the b core. It 
can be shown to be 



V A(cos Ob), (87) 
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where 9h is the angle between r/,i and r/,2. 

After doing the same simplifications as in the light-ion case, we find 

<¥^-^(k,e)|y(rM)|^;"\U,z)) = — i— ^ V(27/ + mil + 1) 



X exp[/(crf + a\)\Pi{cos e)(2l + 1) 

r (88) 

X I draAd/3dydri2drhidrb2raAsin/3ri2rh\rb2 

X P/^cos 0a)A(cos 6b)uifjj.(rAi)ui^jf(rA2)V(rbi) 
X Ui.j.(rhi)ui,jXrb2)MrbB)gi(raA)Pi(cos 6if)/rbB- 

A.l. Coordinates used in the calculations 

We must find the expression of the variables appearing in the integral as func- 
tions of the integration variables rip, r2p, ru, rAa,P, J as indicated in Fig. [19] (re- 
member that Yau = rAa z)- Tau being the center of mass coordinate, we have 

rAa = — (fAl + rA2 + mbYAb) , (89) 



and 



TAa = —(.rAi+ rA2 + mbYAb) , (90) 
ma 

di = — (mfor/,2 -(m^ + l)ri2), (91) 
rtia 

di = — Jinib + lyl + nibinib + l)rl - nibvl^, (92) 
nia ' 

d2 = — ^jinib + lyl^ + nibinib + ly]^ - mbrj^, (93) 



and 



nia niB + mh 

YBb = rAa (fAl + rA2). (94) 

m-b tnbniB 



The angle a between di and ri2 is such that 

cos(ff) = ' . \ (95) 
2rud\ 

The complete determination of rAi,rA2,ri2 can be obtained by writing their ex- 
pression in a simple configuration, in which the triangle lies in the jcz-plane with 
di pointing along the positive z-direction, and r^a = 0. Then, a first rotation K^iy) 
of an angle y around the z-axis, a second rotation Kyifi) of an angle /3 around 
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the i/-axis, and a translation along will bring the vectors to the most general 
configuration. In other words, 

ri2 =5^,W,(7)r'i2, (96) 

with 



1*12 = ni 






di 

sin(Q') 


- cos(a) 



(97) 



(98) 



and the rotation matrixes are 
and 



cos(J3) sm(J3) 

1 
-sm(J3) cos(J3) 

cos(y) - sin(7) 
sin(y) cos(7) 
1 



Vai = 

rAa + di cosQ3) 

d\ sin(/?) + ri2 cos(J3) cos(y) sin(Q') - ru sin(yS) cos(Qr) 
Tai = ru sin(7) sin(Q') 

rAa + dy cos(jS) - sin(/S) cos(7) sin(Qf) - cos(Qr) cosifi) 

We we also find 



and 



We easily obtain 



and 



n\ = —{yai + {mb + l)rAi - m^r^a), 
nib 

rb2 = — (r^i + {mb + l)rA2 - maVAa). 

Mb 



COS 9a = 



'a\ ^ ' A2 ' 12 

2?"ai'"a2 



(99) 

(100) 

(101) 

(102) 

(103) 
(104) 

(105) 
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Figure 19: Coordinate system used for the simultaneous amplitude. The angle a (=0 in 
the figure) coiTesponds to the orientation about the z-axis. The origin is locat ed in the 
center of mass of A, and R is the distance to the center of mass of the triton (from Bay man 

(Eizll))- 



B. Successive and non-orthogonal contributions 



We write the successive transition amplitude as (see lBayman and ChenI (|1982i) ) 



KM 



X [(/r^/(r^i, (ri)i^^/(r^2, (^2)]°* vim) [iA^>(rA2, o-2)i/^Hrbi,cri)]l (107) 

X G(rc,, r^,) [^^/(r:,^^ t^2)'A''(rM > t^'i)]^* ^^ci) 
x[^^'(r;i,(r;),A^'(r;2,cr^)]°y^^(r;j. 

We expand the Green function and distorted waves in a basis of angular momen- 
tum eigenstates: 



/~'>*O^Bb, VBb) = J] H+l-^i-'e-'^fFj(rBh) [Y'(hb)Y'(kBb)l , (108) 



kBb^Bb 



/^\<a) = 2 i' V27TT-i^e''^F;(r;j [Y'(r'^JY'(U]l , (109) 
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where we have taken into account that = z- The Green function in the inter- 
mediate channel can be written as 



G(rcc,rc,) = i }, V24 + 1 7 ; [Y'ircc)Y'(rc,)\ . (110) 



Finally 



= X"TV Z /-K^/'-'V(2/+i)(2/. + i)(2/+i) 



(TlO"2 ^ 

X \Y'ir^.)y'iKl [i"'(fc,)i"'(?y]° ^Jii^l^tl 

X — r ' (fAi, (^i)iA^^ (rA2, 0-2) I 

X [(A^'(r;i, (r;)(A^'(r;2, cr',)]l ^ [<A^>(rA2, (r2)«A^'(rM, (Ti)]^ 



KM 



K* 
M 



(111) 
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Let us now perform the integration over Vai- If we do the approximation of con- 
sidering vai « rci, rAi becomes independent of Tai, so 



^ j dvAi [iA^'^(rAi, o-i)iA^'^(rA2, 0-2)]°* [iA^^(rA2, 0-2)1//^' (rti,o-i)f^ 

= ^(-1)1/2— +i/2-<r, r drA2[>J^^'(rAu-(r,)iff^f(rA2,-cr2)]l[>Jf^^(rA2,^ 

0-1,0-2 

= - r jrA2[<A''^(rA2,-(r2)(A^'^(rAi,-cri)]°[<A^'^(rA2,cT2)<A^'(^^ 

0-1,0-2 ^ 
= -{UfjMjfjdK\UfjMjfjdK)K + 

(T] ,(T2 

X JrA2 [lA^'^ (rA2, -cr2)i/r^>(rA2, 0-2)]° [lA^'^ (rAi , -(ri)(A^'(rM , (tOJ^ 

0"! 



^((/ji)y/Z4),J(Z//,0jf(ii)o)j,[y'^(rAi)F''(^^ 



27/ 

(112) 



where we have evaluated the 9j symbol 



{Ufjf)oUfji)K\Ufjf)oUfji)K)K = T \ 1 • (113) 

^i/ + 1 
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We proceed in a similar way to evaluate the integral over r^j, considering now 
r^2 ~ so r^2 becomes independent of r^^ : 



M 



K 

-M 



= -(-1)''^''(0>j;)^Oj;)oIO>j/)^(jj-,)o)^(- V2j.- + 1) 



M 

K* 



(114) 

Putting all together 
vv ^picci^nfi 2 



T 

^ 2NT 



n kAakBhkcc yjilji + l)(2jj + I) fj;^ 
X 2 ^(2/,+ 1)(2/+ 1)(2/+ 1)?'"' 



X J" £rccd^ri,id^r'c^d^r'^2v(rbi)v(r'^2)^if(rAi)uiXr^^^ 



r'A/Bbrccr'c, 

X [Y\rBb)Y\kBb)'\^ [YKr'^a)YKha)^^ [Y\rcc)Y\r'cS\^ . 

(115) 

We can write 

[Y\hh)Y\ksh)'\^ [y'(r;jy'(^^,)]|| = 

((ll)o(llMllUll)o)o [Y\rBb)Y'(r'^j]l [Y'(kBt)Y'(kAa)]l (116) 

Taking into account that 

\Y'(kBb)Y'ikAa)t = ^Yi,ikBh)i', (1 17) 
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and 



{(llUlc IcMIIcM11c)k)o {[Y'irBh)YHrcc)f [Y'(r'^a)YHrcc)f 



2K+1 



X 



'(2/+ l)(2/,+ 1) 



'(2/+ l)(2/,+ 1) 
X J] [Y'irs,)YHrcc)t [Y^r^JY^rm 



M' 



(118) 



It is important to note that the integrals 

J dfccdhi [Y\hh)Y'^{rce)\^^ [Y'f{hi)Y\hi)X^ , d 19) 

and 

J Jr^,<2 [Y\r'^a)YHr'cc)]l V'' ^^'Ai)Y\n,)]l , (120) 
over the angular variables do not depend on M. Let us see why with (II 191) . 

[Y'ihi,)YHrcc)f^ [Y'firAi)Y'^ihi)]l = (-if''' [Y'irBh)YHrcc)f_^ 

x[Y'f(rAi)Y\hi)]l = {-lf'''Y,^KKM - M\J 0) ^^^^^ 

X [[Y'{rBt)YHrcc)f [Y'f(rAi)Y'-{hi)f}'^ ■ 
After integration, only the term 

(-if-^'iK KM - M\0 0) i^[fCrBb)Y'^(rcc)f [y'^ irAi)Y'- (hi)]"]'^ = ■ 
-^l^= {[Y'(rBh)YHrcc)f [Y'^{rAr)Y'-{hi)f]'^ 



(122) 
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corresponding to 7 = survives, which is indeed independent of M. We can thus 
omit the sum over M and multiply by (2K +1), obtaining 

n^kAakshkcc yj(2ji + \)(2jf+ 1) 

X J](2K + mA)jr(liM¥M\i)ofK (123) 



K 



Ha-'+cr') 

e • t 



with 



Sxik = J d^rccd^rbiv(rbi)u,^.(rAi)ui.(rbi) 



and 



rcc rBb (124) 



Jrcr fixed 



'rccfixed ^Aa ^Cc 

X [t^{r\,)Y\r',,)f* [YHr'c.)y\rA^)]l ■ 

(125) 

The integrand in (I124|) can easily seen to be independent of M, so we can sum 
over M and divide by {2K + 1), to get the integrand 

1 , X , X , SK,uXrcc)Fi{rBb) 

—r—rv{rbOu,{rAi)uiXrbi) 

2K+1 rcc rBb 

M 

This integrand is rotationally invariant (it is proportional to a spherical 
tensor with L = 0, M = 0), so we can just evaluate it in the "standard" configura- 
tion in which Tqc is directed along the z-axis and multiply by 8;:^ (see Fig. I^Dl) . 
obtaining the final expression iox S k,i,Ic'- 

Skjj. - ^^^^ I ^ 

X J rl^drccrl^drbi sin 6 dev(rbi)ui^irAi)ui,(rbi) 
SKiiXfCc) F,{rBb) 

X 

^Cc fBb 

X _^</, / M\K M) [Y'f{rAi)Y\e + n, 0)]^ y^(rB,). 

M 
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Similarly, we have 



SK,uXrcc)= ^^^^ r 



FiKJfiSkcc,r.)PiSkcc,r>) ^^^^^ 



X ^<Z, 01M\KM) [Y'fir'^,)Y'-Crl,)f^ Y'^ir'Aa)- 



{h0lM\KM)[rfir^2)Y-(,r,^y''* 

M 

We obtain the final expression 



^yy _l024iuccn''^i 



2NT 



K^kAakBbkcc -^{Iji + l)(2j/ + 1) 

E ^^((^/5)i/(^4)i.l(¥0^(H)o& (129) 

K 



X 



with 



SKj.k= J rl^drccrlxdrhi smedev{rM)u,f{rci)u^{rbx) 



and 



X '-^^^^^^'^^^P'^'-^^ (130) 
X YjHc / M|if M) [Y'f{rc^)Y\e + n, 0)]^ Y%{hb\ 

M 



SKjjXfCc) = J r'c, dr'c, r}^ dr\^ sin 0' d& v{r[^ui^.{r\^ui.^{r^^ 

' Aa ' Cc 



X Y,ilc / M\K M) [Y'f{r^^)Y\r^^)f^ YUr^,). 

M 



If we are dealing with light ion reactions with spin-orbit terms in the optical po- 
tentials, the amplitudes for diff"erent values of the total angular momentum j of the 
relative motion of the nuclei are combined in the same way as for the simultaneous 
amplitude. 
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B.l. Non-orthogonality term 

We write the nqn-ort hogonality contribution to the transition amplitude (see 



Bayman and ChenI (Il982h '): 



T2nt='^Yj r d'rccd'r,,£rj,2d'rld'r\^x''^*(\ish,rBh) 



0-10-2 

o-'jO-^ 
KM 

X 



(132) 



[(A-'''(rAi, cri)iA-'''(rA2, cr2)f* v(rhi) [(A-'''(rA2, cr2)iA-'''(rM, cri)]|^ 

This expression is equivalent to (11071) if we make the replacement 

^G(rc„ r'cMKi) ^ 3(rcc - r'c,). (133) 

Looking at the partial-wave expansions of G(rcc, r^^) and 5(rcc -r^c)' ^^^^ 
we can use the above expressions for the successive transfer with the replacement 

i-r2 ; K^c2) ^ ^(rcc - rcc)- (134) 



We thus have 



with 



'""^ ^A^^B/, V(2j- + I)(2j> + 1) 

^ 2^^^^2)-'>*'^'2^^'''(^/^'^^*^H)o& (135) 

SKj.h= J rl^drccrlidrbi sin 6 dO v(rhi)uif.(rci)ui.(rhi) 

x^£^^ (136) 

fCc fBb 

X _^</, / M\K M) [Y'f{rci)Y\e + n, 0)]^ Y'^(rBbX 



M 



and 



SKjjXrcc) =rcc \ dr'j^^rli sinO' dO' Ui^.(r^2)Mrc2)-^'^"^ 



Aa 

K* 



X ;^</, / M\K M) [YHf\2)Y\r[2)\* yl,(rl) 



(137) 



M 



B.2 Coordinates used in the calculation of successive transfer 



61 



B.2. Coordinates used in the calculation of successive transfer 

In the standard configuration in which the integrals (11301) and (11311) are to be 
evaluated (see Fig. [2D1) , we have 



Now, 



rcc = rcc^, r^i = rii(-cos6iz-sin0x). (138) 



_ _ 

Tci - Tcc + Tel - TCc "I T^b\ 

mi, \ ^ f^b 

rcc -rh\ cose z -rbi sinfe, 

Mb + I I nib + I 



and 



But 



and 



(139) 



rsb = Tfic + rcb = — —rci + rcb- (140) 
niB 

rcb = rcc + Vcb = rcc ^rM, (141) 

mb + I 

so, substituting in (11401) we get 

/ms-l lUb + mB \ Mb + mB . 

YBb = rcc + — ; —rbi cos 61 z + — ; — sin 9x. (142) 

\ Mb niBimh + 1) / mB{mb + 1) 

The primed variables are arranged in a similar fashion, 

r'c, = r'c,z, = r'^i- cos ffz- sin ff x). (143) 

And we get 

^c2 = -fee - zr—TfAi COS z - — — rr^2 sm e x, (144) 
\ + 1 / + 1 



^cc ; 7T''a2Cos6' z —r^2^m.ex. (145) 

m„ ma{mA + IK / '^.X^a + 1K' 
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Figure 20: (a) "Standard" configuration, with r^c along the z-axis, r^i in the x - z plane 
and {ri,i)x < 0. ( b) The result of applying a rigid-body rotation {a,p,y) to the configura- 
tion of (a) (from Ib ay man and ChenI (Il982[) ). 
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